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Abstract 

For the first time a method for reahzing macroscopic quantum optical soh- 

tons is presented. Simultaneous photon-number and momentum squeezing is 

predicted using soliton propagation in an interferometer. Extraction of soliton 

pulses closer to true quantum solitons than their coherent counterparts from 

mode-locked lasers is possible. Moreover, it is a general method of reducing 

photon-number fluctuations below the shot-noise level for non-soliton pulses 

as well. It is anticipated that similar reductions in particle fluctuations could 

occur for other forms of interfering bosonic fields whenever self-interaction 

nonlinearities exist, for example, interacting ultracold atoms. 
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Quantum solitons are energy eigenstates or photonic number states for optical systems. 
The particle number fluctuations of these quantum states are zero and this property is 
preserved by the nonlinear system. A quantum soliton is a fundamental state in nature 
capable of allowing information bits to propagate arbitrarily long distances in the presence 
of dispersion and nonlinearity. Even approximate forms of such an idealised quantum object 
have been difficult to observe. 

Mode-locked lasers do not produce quantum solitons but instead have a Poissonian dis- 
tribution of number states Solitons excited by such laser pulses have initial fluctuations 
in the four soliton parameters of photon number, momentum, position and phase. Neither 
are they minimum uncertainty states. Position fluctuations grow quadratically for a freely 
propagating fundamental soliton due to the linear dispersion acting on the initial momen- 
tum fluctuations. Similarly phase noise grows due to the initial fluctuations in its conjugate 
variable photon number acted upon by the nonlinearity. Typically initial fluctuations in 
position and phase become insigniflcant compared to the increase due to quantum diffusion. 

In this Letter quantum soliton generation using an approach based on interference of op- 
tical flelds is proposed. The idea that one can manipulate the internal quantum noise struc- 
ture of propagating initially coherent solitons to produce more than lldB photon number 
squeezing using an asymmetric Sagnac loop was flrst put forward by the author recently . 
That followed surprising results on interference of solitonic flelds in nonlinear optical loop 
mirrors to produce more than 15dB excess noise 0. Here it is shown that in the easily 
accessible regime of macroscopic photon numbers of order 10^ — 10^ typical in picosecond 
and subpicosecond laser experiments, it is possible to produce soliton pulses with greater 
than lOdB reduction in photon number fluctuations. In addition, momentum fluctuations 
can be reduced to more than 6dB below that of a coherent state pulse. Initial experiments 
have recently observed 3.9dB (6.0dB inferred from 79% detection efficiency) ^ (5.7 dB 
has independently been observed) photon-number squeezing at room-temperature using the 
method disclosed earlier and discussed here in more detail. 

Any scheme which involves interference of bosonic flelds (different in either direction 
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or polarization) and where at least one has undergone evolution according to an effective 
(1 + 1)D nonlinear Schrodinger equation could exhibit similar behaviour. This includes 
waveguided atomic solitons [^], spatial optical solitons for nonlinearities that depend on the 
particle density flux and optical pulses in cascaded quadratic media or phase-mismatched 
X*^^^ simulton interactions |^. Results are presented for an ideal Mach-Zehnder interfer- 
ometer since the two counterpropagating fields of a Sagnac fiber loop are assumed to have 
propagated in two independent waveguides. It is important that the system is described us- 
ing a dispersive quantum field theory in order to contain the essential physics and to make 
predictions for current optical pulse experiments which use pulse durations corresponding 
to 70fs < to < 2ps. For optical systems the nonlinearity and dispersion can be turned on/off 
in one arm easily. For systems containing real massive interacting particles the loop results 
are more appropriate. Optical experiments were performed using Sagnac loops to reduce 
low-frequency noise. Importantly, the scheme is not critically sensitive to exact coupling 
ratios, powers or fibre lengths although optimization is required to reach beyond the lOdB 
limit. The interference of optical fields has several distinct advantages over direct spectral 
filtering ^ including greater noise reduction, smooth output pulse envelopes and simulta- 
neous photon-number /momentum squeezing. The model investigated to demonstrate the 
idea - the nonlinear Schrodinger equation - is ubiquitious in the physics of dispersive self- 
interacting fields and describes perturbations in a wide variety of physical systems. It is 
anticipated that particle number fiuctuations in more general systems such as described by 
a quantum Ginzburg-Landau equation could exhibit similar behaviour to the idealized soli- 
tons discussed in this Letter. There are strong reasons to believe this is the case especially 
for weak damping which occurs in optical fibers in the 1.5/im regime. 

The nonlinear Schrodinger equation has been used in various forms for the study of Bose- 
Einstein condensation 0, waveguided ultracold atoms and propagating coherent quan- 
tum optical solitons p0|-p!3| where the quantum nonlinear Schrodinger equation (QNLSE) 
governs dynamics of the photon fiux amplitude. The Raman-modified stochastic nonlinear 
Schrodinger equation for the normalised photon fiux fields {0(C) t), <PKC: ''")} iii the positive- 
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P representation is given by 
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where length and time variables (C, r) in the comoving frame at speed uj' (group- velocity at 
the carrier frequency) in the laboratory frame (x, t) are r = (t — x/uj')/tQ,( = x/xq,Xq = 
to/|/c"|- For this equation and its Hermitean conjugate for (f)\ the characteristic time scale 
to will be chosen to be the pulse width and the soliton period is 7r/2 times longer than the 
dispersion length xq determined by to and the second-order dispersion k". The quantum noise 
from the electronic nonlinearity Fg is a real delta-correlated Gaussian noise with variance 
given by the product of the electronic fraction / (ideal QNLSE has / = 1, /i(r) = F^, = 
0) and inverse photon number scale 1/n = xto / \k" \uj''^ ■ Silica fiber Raman gain has a 
peak near 13 THz and / = 0.81 for Raman inhomogeneous model parameters used in 
evaluating the response function /i(r) and noise F^ corresponding to the Raman gain curve 
in Reference [Q. Quantum field propagation is performed numerically using techniques 
discussed elsewhere |15|. All simulations without Raman used n = 10^ with averaging over 
10^ trajectories for the positive-P representation. Error bars in the plots represent the 
estimated combined sampling and step-size error. 

To illustrate the interference at the output beamsplitter of unequal amplitude solitonic 
fields consider the interference term of the transmitted photon spectral flux proportional to 



The interference term above is the same as obtained from a generalised quadrature-phase 
operator measurement in homodyne detection where the local oscillator (weak field) can 



a state preparation device and later, combined with the weak field and photon detector, as 
a measurement apparatus. It is easy to show that placing a phase-shifter in the weak field 
arm can be used to switch from sub-Poissonian to super-Poissonian statistics in analagy 
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have nonclassical statistics [|T^. One can consider the beamsplitter of a loop to act first as 



with changing a quadrature-phase measurement from the squeezed quadrature to the anti- 
squeezed quadrature. Since the two fields 0i and 02 in each arm of the interferometer 
typically has a different spectral distribution due to the different initial amplitudes, there 
exists a spectral filtering mechanism. Importantly, it has been shown previously that the 
internal quantum noise structure for the two fields will be quite different P,p!7|. These 
two effects combined will alter the quantum statistics of the resultant field. In addition, 
the interferometer is obviously sensitive to pulse chirp, a frequency dependent phase-shift, 
induced by group-velocity dispersion which is important for fiber lengths of a few soliton 
periods discussed here. These aspects are not present in any single-mode description as used 



in Reference where single polarization and nonlinear polarization rotation interferometric 
configurations were considered with one arm as free space. This paper goes beyond the one 
arm configuration and shows that noise reduction is possible even for solitonic fields in both 
arms. 

While a heuristic analysis based on the energy input-output curve for Sagnac loops can 
sometimes estimate whether the output photon-number noise is expected to be below shot 
noise, it is generally not reliable. It is suggested that an appropriate method for large photon 
numbers which does not require the more sophisticated techniques used with the positive- 
P representation is to use the truncated Wigner representation |]l9l which for the ideal 



QNLSE with coherent state inputs involves solving only the classical equations argmented 



with Gaussian noise on the initial conditions [15|. The classical nonlinear phase shift and 



dispersion in each arm of the loop, which has long been known to support effective soliton 
switching pO[ and reduction of dispersive waves [^], leads to a characteristic input-output 



curve for the highly asymmetric (90:10) and near balanced (60:40) case. The transmitted 
pulse photon-number (scaled by n) versus the input amplitude with 0(0, r) = A^sech(r) 
inputs is given in Fig. |1| for the case 90:10 at ( = tt, 2tt and 60:40 at ( = 27t. A propagation 
distance longer than C = vr for the 60:40 case was chosen so that input-output curves contain 
at least one turning point. One can see that after propagating 2 soliton periods for the 90:10 
case, the input-output curve's slope is positive leading to excess noise at the output for all 



inputs except near the turning points N = 1.35, 1.5, 1.85 where a saturation effect might be 
expected. The double dip structure in Fig. 2, which has been observed experimentally [Q, 
corresponds closely to these classical turning points. Slightly beyond the classical turning 
points a negative slope represents a region stable against changes in the input state. In 
the 60:40 case, the input-output curve in Fig. |^ suggests that significant squeezing would 
occur near = 1.62 where the slope changes sign. The quantum- field theoretic results 
discussed next do not predict noise reduction for N = 1.62, ( = 27t — in sharp contrast 
to the simplified classical input-output picture the quantum theory actually predicts 15 dB 
excess noise in this case which the truncated Wigner theory also predicts. A simple 
explanation for this disagreement lies in the difference between a loop and an interferometer 
with only one nonlinear pathway. In the former situation quantum theory allows quantum 
phase diffusion to develop in both arms independently whereas the classical argument would 
have input noise appear reduced in strength in both arms. Clearly this is exacerbated in the 
60:40 loop compared to a 90:10 loop. The heuristic argument is expected to disagree with 
quantum theory for a 90:10 loop but with a much reduced error compared to the 60:40 case. 
Quantum field-theoretic results for the 90:10 case will now be discussed in more detail. 

For a fiber loop wih a beamsplitter transmission of 90% large photon-number squeezing in 
the transmitted field (i.e., the "dark" port for a 50:50 beamsplitter arrangement) is predicted. 
Large noise reduction occurs over a wide variation of coupling ratios as well. Although 
optimal parameters for the largest reduction in photon-number fluctuations are still under 
investigation, it is predicted that 11 ± 1 dB squeezing at C = 3 is possible for coherent 
A^sech(r) input pulses in the absence of Raman noise for A^ = 1.5. Obviously smaller energy 
pulses can be used in combination with longer propagation distances. Significant reduction 
in the momentum fluctuations, up to 6 dB below that for a coherent pulse, also occurs but the 
latter result is preliminary. In the A^ = 1.5 case the lower pulse energy arm of the Sagnac 
loop is only dispersive radiation while the other arm contains a solitonic field (emergent 
soliton plus dispersive radiation) whose internal noise structure was recently described by 
T% . As discussed earlier the nonlinearity in the lower energy arm is not required to 
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observe noise reduction but must be included to correctly predict the noise levels expected 
from a loop. 

The predicted variation of the photon-number squeezing in dB versus the input pulse 
amplitude parameter with 0(0, r) = A^sech(r) coherent inputs after propagating in a 
Sagnac loop of length C = using the ideal QNLSE is given in Fig. 0. The influence of 
the Raman effect for to = O.lps at room temperature (shown in Fig. 3) was determined to 
be not large in this case even though the input pulse to the fiber had more than twice the 
energy of a fundamental soliton in the case of = 1.5. The photon- number fluctuations 
are however still increased by the Raman effect to 8.3 ± 0.4dB below shot-noise at C = vr 
for = 1.5 including O.ldB/km losses. There is a strong similarity between variation 
with energy for a fixed distance and variation with distance for a fixed input energy since 
the nonlinearity allows soliton pulses of larger (smaller) energy to experience similar effects 
over shorter (longer) distances. Both in Fig. 2 and Fig. 3 the loop output exhibits excess 
noise, which was also predicted |jl2| and observed for direct spectral filtering, and is 



clearly a general feature in nonlinear systems. After A^ = 1.5, the low energy pulse switches 
from dispersive radiation to a solitonic field with an internal noise structure whose spectral 
correlations can lead to excess noise or squeezing depending on the bandwidth of any filtering 
mechanism |T^J^. The noise reduction predicted for a one soliton period fiber is also given 



to demonstrate the possibility of significant squeezing in short fibers provided that the pulse 
launched into the fundamental propagating mode of the fiber matches the initial conditions 
assumed here. The variation for = 7r/2 appears similar to ( = 7C but with higher energy 
required as expected for the nonlinear fiber. For comparison a much longer propagation 
distance of 16 soliton periods using A^^ = 10/9 is given in Fig. 3 so that a fundamental 
soliton propagates in one arm and the weak pulse experiences free propagation (non-loop 
case). While solitons are not necessary the squeezing predicted for pulses in the normal 
dispersion regime is less and at higher energy than soliton pulses as expected. In this case 
the pulse quickly temporally broadens from the group velocity dispersion while the pulse 
spectrum broadens initially and then reaches an equilibrium even for A^ > 1. This is in 



contrast to the breathing in the anomalous regime. In Fig. 3 the case = 3 is given 
without Raman and reaches 6.4 ± O.ldB below shot noise. 

In summary, photon number fluctuations for the two coupling regimes of the Sagnac 
loop - slightly asymmetric and highly asymmetric - have quite different characteristics. 
The nonlinear optical loop mirror (slight asymmetry) usually increases the photon-number 
fluctuations above the shot-noise level for coherent inputs P]. Sub-shot noise statistics is 
possible however the 60:40 loop is much more sensitive and without precise control would 
almost certainly produce excess noise at the output. On the other hand, photon-number 
noise may be significantly reduced below the shot-noise level easily for highly asymmetric 
coupling. We have shown for the first time that large photon-number squeezing of solitonic 
fields can be produced by a Sagnac loop over a significant range of input energies and the 



squeezing is larger than predicted for spectrally filtered optical fiber solitons |T^. Despite 
this similarities in the output photon statistics for the two approaches exist due to the use 
of the same nonlinear propagation to produce the internal quantum correlations of a quan- 
tum soliton not present in the initial coherent state. In essence a nonlinear interferometer 
is capable of producing optical pulses closer to true quantum solitons than their coherent 
counterparts from mode-locked lasers. The application of these ideas to controlling fluctu- 
ations in atomic wavepackets is particularly interesting given recent experimental progress 
with Na atom coherent condensates with densities ~ 10^^ cm~^ using optical traps p2 
This is several orders of magnitude larger than estimated to observe atomic sohtons H. 
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FIGURES 

FIG. 1. Transmitted pulse photon-number scaled in units of n versus the input amplitude N 

with 0(0, r) = A^scch(r) inputs. The input-output curve is sliown for ^ = vr, 27r using a 90:10 loop 
and also for ( = 2Tr for a 60:40 loop. 

FIG. 2. Photon-number variance (dB) at ( = 7^/2, ( = tt for 90:10 Sagnac loop with 
0(0, r) = A?^sech(r) coherent inputs using the ideal QNLSE. 

FIG. 3. Photon-number variance (dB) versus propagation distance for a coherent 
0(0, r) = A'^sech(r) pulse inputs into a 90:10 interferometer. The horizontal axis is rescaled for 
the case N'^ = 10/9 (solid line) where the lower energy pulse undergoes free propagation with no 
Raman. Sagnac loop results are given for N = 1.5(dashed line) without Raman and with Ra- 
man using to = O.lps (n = 10^) at phonon temperature of 300K. An example of a non-solitonic 
input, A;" > 0, is shown ior N = 3 without Raman. A''^ = 10/9 and N = 3 cases used truncated 
Wigner theory averaged over 5000 trajectories. 
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